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Abstract: Recent CMS data significantly extend the exclusion limits for supersymmetry. 
We examine the impact of such data on global fits of the constrained minimal super- 
symmetric standard model (CMSSM) to indirect and cosmological data. By simulating 
supersymmetric signal events at the LHC, we construct a likelihood map for the recent 
CMS data, validating it against the exclusion region calculated by the experiment itself. 
A previous CMSSM global fit is then re-weighted by our likelihood map. The CMS results 
nibble away at the high fit probability density region, transforming probability distribu- 
tions for the scalar and gluino masses. The CMS search has a non-trivial effect on tan /3 
due to correlations between the parameters implied by the fits to indirect data. 
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1 Introduction 

Many recent efforts have examined whether various simple supersymmetric models are com- 
patible with the anomalous magnetic moment of the muon, the dark matter relic density, 
direct searches for supersymmetric particles and Higgs bosons, and electroweak observables 
simultaneously [1-18]. Global fits take into account variations with respect to all of the 
parameters of the model that have a significant impact on the observables, including Stan- 
dard Model (SM) parameters. Various algorithmic tools have now been developed to allow 
such a sampling of a multi-dimensional parameter space, which may be multi-modal [19- 
21]. However, current global fits of indirect and cosmological data to the CMSSM are 
not robust. This is shown by the large prior dependence of the Bayesian fits [2, 5, 7, 10]. 
Prequentist fits to edge measurements from hypothetical LHC SUSY edge measurements 
showed an incorrect confidence level (C.L.) coverage of frequentist fits when the C.L.s are 
calculated by assuming a distribution [22]. There are no published coverage studies of 
global SUSY fits and, since they are expected to be less robust than fits to an LHC SUSY 
signal, a coverage study of the frequentist fits is necessary and long overdue. A fit to a 
large volume string model with only two free parameters additional to the SM (the ratio of 
the Higgs vacuum expectation values, tan /3, and an overall supersymmetry breaking mass 
scale) did display approximate prior independence [23]. On the other hand, a fit to a model 
with three parameters additional to the SM (minimal anomaly mediated supersymmetry 
breaking) showed significant prior dependence [24]. Fits to models with more than three 
additional parameters have also (so far) shown a lack of robustness [11, 13, 24, 25]. 

The LHC general purpose experiments ATLAS and CMS actively searched for super- 
symmetric states in 2010, collecting some ~ 35pb~^ of recorded collisions each at a center 
of mass energy = 7 TeV. The first dedicated analysis of this data recently appeared 
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from CMS, examining multi-jet states accompanied by missing transverse momentum [26]. 
The data were statistically compatible with the SM predictions, allowing the experiment to 
place constraints on supersymmetry. In its publication, CMS also analyzed the CMSSM, 
showing that the equivalent previous best exclusion from experiments at LEP and the Teva- 
tron are significantly surpassed. In 2011, some thirty times more luminosity is expected 
to be collected by the LHC experiments than in 2010 and, in the absence of a signal, the 
search reach will doubtless be extended significantly. 

Here, we wish to examine the extent to which the first publicly available search from 
CMS impacts on the previous global supersymmetric fits. Despite the fact that such fits 
are not yet robust and therefore not definitive (and are likely to remain so until direct 
supersymmetric searches yield a significant signal), such an exercise should give interest- 
ing information on whether the good-fit region is being covered yet, and the prospects for 
searches in the near future. The effect of the CMS direct SUSY search on the fits may 
be qualitatively examined. We demonstrate how one can perform such an analysis in a 
reasonable amount of computational time, and hopefully set precedents for good practice, 
such as validation against CMS's more sophisticated event simulation. The CMS search 
has already been used to examine how much of the CMSSM parameter space is allowed 
beyond the 95% exclusion contour: [27], around 1% according to the author's metric related 
to the naturalness of electroweak symmetry breaking. We shall examine the effect of the 
search results on a previous Bayesian fit, because of the computational ease with which 
the impact of the direct search can be quantified (as opposed to the frequentist methods, 
where an examination of coverage would be necessary to set confidence levels, requiring 
many such fits). We pick the CMSSM as our model partly because CMS analyzed it ex- 
plicitly in their publication. Thus, it is possible to check our evaluation of the CMS data, 
which includes various approximations and simplifications, against their more sophisticated 
treatment. We use the 95% confidence level exclusion contour in a two dimensional param- 
eter plane that CMS provided in their paper. The CMSSM is a familiar model that many 
works have analyzed and it provides a well-defined playground for examining supersym- 
metric phenomenology, with only a few free parameters additional to the SM: the universal 
scalar mass mo, the universal gaugino mass mi/2i the supersymmetric (SUSY) breaking 
scalar trilinear coupling Aq, tan/3 and the sign of fi, a parameter in the Higgs potential. 
Although the CMSSM is very specific and may well not represent the correct pattern of 
supersymmetry breaking, aspects of its phenomenology are similar to the class of models 
which are effectively perturbations around the CMSSM assumptions. 

Early ATLAS jets plus missing transverse momentum data presented at Summer con- 
ferences in 2010 have already been used [28] to place bounds on simplified models containing 
squarks and gluinos which decay to jets and neutralinos. This analysis was based only on 
70 nb~^ of data. This data is not used here because it is only a small integrated luminosity, 
with a small effect compared to the data we use. We shall see that the predicted CMS 
SUSY search signal is approximately independent of all parameters except mo and mi/2, 
which leads to significant simplifications in the incorporation of its results into the fits. 

Our paper proceeds as follows: in section 2, we detail the CMS SUSY search and 
results, showing how we approximately reproduce their analysis, resulting in a likelihood 
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map on the CMSSM parameter space. We then go on to apply the approximate hkeUhood 
map to global CMSSM fits in section 3, showing how the fits change. We conclude in 
section 4. We calculate the accuracy of our approximation that the CMS SUSY signal is 
approximately independent of all parameters except niQ and mi/2 in Appendix A. 

2 CMS ar Likelihood Map 
2.1 The ar Search 

In 35pb~^ of pp collisions at ^/s = 7 TeV at the CERN LHC, CMS examined events 
with significant transverse momentum |pt| in the jets ji for multi-jet events, i.e. Ht = 
Si=i* IPtI '^^^ GeV. pj, = {px, Py) is the jet momentum transverse to the beam. 
Isolated lepton and photon vetoes were also applied. The final event selection relies on a 
variable ot, that is designed to discriminate effectively against QCD multi-jet production, 
where one of the jets' transverse momentum is significantly mis-measured [29], although 
it has been argued that the variable gives no special immunity to initial and final state 
radiation effects [30]. The system of Nj jets is reduced to a system of two jets by combining 
them into pseudo-jets A and B^. The combination chosen is the one that minimizes 

A/?t=EIPtI-EIPtI- (2.1) 



One then calculates 



"T = — , ^ , (2.2) 



where = V (Sfif Pi")' + (EI =1 PyY is defined at the jet-level. Note that often, a 
different definition of q;t is given which is identical in the idealized case of massless di-jet 
events, but which differs from the result of Eq. 2.2 for the case of more than two jets. The 
cut qt > 0.55 is used in CMS's analysis. 

The jets are defined by the anti-A;^ algorithm with a size parameter of 0.5, a minimum 
IptI > 50 GeV and a bound on the pseudo-rapidity 77 of |r/| < 3. The highest |pt| jet 
is required to have \r]\ < 2.5 and the transverse momentum of the two leading jets must 
exceed 100 GeV. T^rp[cal) is defined over all calorimetric deposits of the event, and events 
with i?niiss ~ -^t/Ii^t('^"OI > 1-25 are rejected to protect against multiple jets failing the 
IptI > 50 GeV requirement. 

Additional experimental cuts to remove events where significant energy losses from jets 
going into un-instrumented regions were also placed. After all of these cuts, o = 13 events 
on a predicted SM background of 9.3±0.9 were observed. This search is compatible with 
SM predictions at the 2a level, with a ~ la slight excess of events. 



^ We note here that the difference between transverse momentum and transverse energy can be significant 
for pseudo-jets, even within an event generator, because they have effective masses. Like the authors of 
Ref. [31], we urge care when defining them. 
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2.2 Simulation of the SUSY signal 

In order to simulate the production of sparticles from LHC pp collisions, we use HERWIG++-2 .4.2 [32] , 
with the default underlying event model switched on. Supersymmetric spectra are gen- 
erated with S0FTSUSY3 . 1 . 7 [33], with decay branching ratios and widths calculated by 
SDECAY [34]. The information about supersymmetric masses, mixings, couplings and decays 
is transferred between the programs by the SUSY Les Houches Accord [35] . f ast j et-2 .4.1 [36] 
is used to define the jets in the anti-Zcr scheme with R = 0.5 in the energy recombination 
scheme. The cuts mimic those of the experiment and are summarized in Eq. 2.3: 

Ht > 350 GeV, \p^ \ > 100 GeV, ar > 0.55, iJj^igg < 1.25, / isolation. (2.3) 

The lepton isolation veto (7 isolation') is implemented as follows^: any events with leptons 
that have less than 10 GeV of \pt\ in a cone of AR = \J (A?;)^ + (A0)^ < 0.3 are vetoed, 
where (\) is azimuthal angle around the beam-line. For a given CMSSM parameter point 
predicting s supersymmetric events passing the cuts, we calculate the likelihood to be 

C = (2.4) 

from the Poisson distribution for the central value of the predicted number of signal s 
plus background events h = 9.3. We neglect the small uncertainty of 0.9 on h: in the next 
section, where we validate our calculation, we shall see that the effect of our approximations 
(including this one) are reasonable. 

Our simulation of the signal is cruder than the CMS analysis: it is only leading order 
in QCD, whereas CMS's includes next-to-leading order factors for the overall cross-section, 
and we have not performed a dedicated detector simulation to convolute the signal kine- 
matic distributions with the detector response. Even had we included these efi^ects, it 
would be important to check our approximations by reproducing CMS's calculated signal 
because, being outside the collaboration, we do not have access to the full detector simu- 
lation. There is no need for us to validate the SM background, since we may use CMS's 
predicted rates for it. 

2.3 Validation of the SUSY signal 

In order to validate our calculation of the signal and see that we are obtaining a reasonably 
accurate result despite our approximations, we first calculate the ot distribution in the 
di-jet and more than two jet channels for the CMSSM model point LMO (mo = 200 GeV, 
mi/2 = 160 GeV, Aq = -400 GeV, tan/3 = 10 and /U > 0). The only cut we place for 
the purposes of this check is Ht > 350 GeV, mimicking CMS. The distributions for the 
'di-jets' and 'three or more jets' channels are shown in Fig. 1 for our simulation and for 



^CMS uses a lepton isolation criterion involving tracks and energy deposits in calorimeter cells. Since 
we do not perform a detector simulation, we do not have access to simulated tracks and so we must deviate 
slightly from CMS's lepton isolation criteria, although our criteria are very similar in effect. CMS' lepton 
isolation criteria are: in AR < 0.3 around the lepton, the lepton is considered isolated if the scalar pr sum 
of all of the tracks in the cone plus the |pt|s of the energy deposits in the calorimeter divided by the \pt\ 
of the lepton candidate is larger than 0.15. 



-4- 



'.Allaiiach.arXiv:! 102.3149 



■a 




0.01 



0.2 0.4 0.6 0.^ 



1 1.2 1.4 



Figure 1. SUSY signal ut distributions for the LMO point and ^/s = 7 TeV in the di-jets channel 
and in the 3 or more jets channel, as displayed by the legend. Solid lines are obtained from Fig. 2 
of Rcf. [26]: the CMS signal simulation including next-to-leading order (NLO) corrections and full 
detector simulation whereas the dashed lines show the results of our simulation and approximations. 
The only cut applied is Ht > 350 GcV. 
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Figure 2. Our approximation to the CMS ar-search CMSSM likelihood map for (a) tan/3 = 3, 
Aq = 0, (b) tan/3 = 30 and Aq = —200 GeV, where the blacked out region at the bottom denotes 
a fi lightest supersymmetric particle. The region below the red (lighter) curve is excluded at 95% 
confidence level (C.L.), and Ax^ is clipped at 10. The CMS 95% C.L. curve is shown as the blue 
(darker) line. 



CMS' signal simulation. The figure verifies that our calculation of the ay distribution is 
compatible with the calculation of the experimental collaboration: the normalization of 
the sub-dominant exclusive two jets sample is slightly different, but the shapes of both 
samples match extremely well. After all of the other cuts are applied, including ot > 0.55, 
the acceptance times efficiency of the SUSY signal selection is 5.0%. 

Next, we perform a scan over CMSSM parameter space to see how closely we can 
reproduce CMS's calculation of the 95% contour. Like CMS, we choose tan/3 = 3 and 
^0 = for this scan. At each point in an 11 by 11 grid, we simulate 10 000 SUSY events. 
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calculating C from Eq. 2.4. We then calculate Ax^ = — 2 ln£/£(max), where £(max) is 
the maximum value of the likelihood over the plane. A^^ is then interpolated between the 
grid points, and shown in Fig. 2a. The background colour density displays Ax^, and is 
clipped at an upper value of 10. The light green solid curve is our predicted 95% C.L. lower 
bound, coming from the Ax^ = 5.99 contour. We compare this with the blue (darker) line, 
the published CMS contour for the full next-to- leading order observed limit. We see that 
despite our approximations, the contours for the expected sensitivity and the 95% C.L. 
limits agree reasonably well. Indeed, next-to-leading order corrections to the production 
cross-section will tend to increase the signal rate, whereas we expect detector effects to 
decrease it. Thus we expect that the two effects cancel to some extent. We conclude that 
we should be able to use the likelihood as calculated here, although we note that we would 
really like more likelihood information from the experiments far away from the 95% level 
contour so that we could validate our simulation more thoroughly. The ~ \a excess of 
observed events manifests as a dark valley of < 1 in Fig- 2. This will have an effect on 
the CMSSM global fit, as we shall see. 

For the ut search, where the signal involves just jets and missing transverse momen- 
tum, we expect the signal rate to be approximately independent of tan/3 and ^o- This is 
because the signal is dominated by the strong cross-sections of squark and gluino produc- 
tion, which do not depend to any significant degree on those parameters. Third generation 
squark masses do depend upon and tan/3 to some extent, but the SUSY production 
cross-sections have a dominant component coming from the two lighter generations of 
squark and the gluino. The decay cascades of the first two generations of squarks and 
gluinos are likewise expected to be quite insensitive to and tan /3, being dominated 
by strong processes until the decay into the lightest neutralino. Being able to model the 
dependence of the CMS a^-search likelihood on mo and mxi2^ while ignoring the effect 
of ^0 and tan/3 leads to a significant simplification when we come to take it into account 
in our global CMSSM fits. We check the assumption of independence of £ with respect 
to Aq and tan/3 with another mo-mi/2 scan at a different parameter point: = —200 
GeV and tan /3 = 30, for illustration. The result is shown in Fig. 2b: the likelihood is 
similar to the one in Fig. 2a, and changing Aq and tan /3 has not had a significant effect. 
Other constraints on parameters, such as those coming from having a charged lightest 
supersymmetric particle (shown by the blacked out region at the bottom of Fig. 2b) do 
display a significant tan /3 and dependence, but these are already taken into account 
within our fits and will not pose a problem. The very bottom of the plot has negative mass 
squared staus, and the signal rate could not be reliably calculated there. This has a slight 
effect on the interpolated 95% C.L. contour, and is responsible for it having an apparent 
(false) turning point in mi/2- We shall quantify the effects of the tan (5 — Aq independence 
assumption in our fits in Appendix A. 

For increasing mi 12-, we see Ax^ reaching a constant in Fig. 2 because there is no 
SUSY signal, since squarks and gluinos become too heavy to be produced. At large rriQ 
and small the SUSY signal is strongly dominated by gluino pair production, where 

the gluinos have three-body decays into squarks. Thus the dependence of C on mo, if it is 
above 1400 GeV, is negligible. We shall therefore model the likelihood as follows: we use 
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s = for mi/2 > 500 GeV. For mo > 1400 GeV and mi/2 < 500 GeV, we use the C value 
given by the mo = 1400 GeV hue on the figure. For mo < 1400 GeV and mi/2 < 500 GeV, 
we interpolate linearly within the grid of likelihoods calculated. 

3 Global CMSSM Fits Including the CMS ar Search 

We shah use a previous global Bayesian fit of the CMSSM from the KISMET collaboration [5] 
to: the relic density of dark matter, the anomalous magnetic moment of the muon, previous 
direct searches for sparticles, the branching ratios BR{b — t- 57), BR{Bs — s- ^/j), My/, 
siv? 6\^j, as well as 95% exclusions from LEP and Tevatron direct search data. In order 
to predict BR{b — )• 57), BR{Bs — ?• fJ,fi), the dark matter relic density and the anomalous 
magnetic moment of the muon, micrOMEGAsl . 3 . 6 [37, 38] was employed. It computes 
BR(b — )• 57) including one-loop SUSY corrections and NLO QCD corrections. BR[Bs — )• 
/i/i) is calculated to one-loop. micrOMEGAs calculates the non Standard Model correction to 
the anomalous magnetic moment of the muon to one loop, and we added the dominant two- 
loop SUSY QED contributions from Ref. [39, 40]. SUSY-POPE [41], which includes SUSY 
contributions up to two loops, was used to calculate My/ and sin^ 6*^. Spectral predictions 
from S0FTSUSY3 .1.7 are used to place bounds from pre-LHC direct searches for sparticles 
and the lightest CP-even higgs boson. S0FTSUSY3 . 1 . 7 calculates the sparticle masses to 
two-loop in the renormalisation group equations, and one- loop in the threshold corrections. 
The dominant two-loop MSSM corrections are added to the Higgs mass computation, and 
a 2 GeV error on the prediction of m^ was assigned. 

The ranges of input parameters considered were: 2 < tan/5 < 62, |^o|/TeV < 4, 
60 < mi/2/GeV < 2000, 60 < mo/GeV < 4000. Variations with respect to the top 
mass mt, the strong coupling constant as{Mz)-, the fine structure constant a{Mz) and the 
bottom quark mass mb{mi,) were all included. Various different prior distributions were 
examined in Ref. [5], but here we shall use the example of priors flat in the parameters 
listed above, except for mo and mi/2, which are flat in their logarithm. Using such log 
priors allows us to illustrate the effects of the ax search more acutely than with purely 
linear priors. Rigorous convergence criteria were satisfied by the fits, which were performed 
by ten Metropolis Markov Chains running simultaneously. For more details on the fits, we 
refer interested readers to Ref. [5]. 

The effect of the data is encoded within a likelihood function L{m) defined on model pa- 
rameters m, and Bayesian inference consists of turning a prior probability distribution 7r(?n) 
into a posterior probability distribution p(m|(i) via Bayes' theorem p{m\d) = L{m)TT{m)/ D. 
D = J dm L(m)7r(m) is a constant over parameter space: the Bayesian evidence. Since 
we are performing parameter estimation, we are never interested in normalizing constants, 
because the inference is performed under the hypothesis that the model we are examining 
is correct and therefore the constants are always normalized such that the total posterior 
probability is 1. We take 2.7 million points from the Markov chain fits. The density of 
sampled points is proportional to the posterior probability density at any point. Each point 
i typically appears sampled several times, and we therefore set its weight pciPi) equal to 
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Figure 3. Global CMSSM fits: (a) excluding the CMS ar search and (b) including the CMS 
ut search likelihood. The posterior probability of each bin is shown as the background colour, 
normalized to the maximum bin probability. The almost vertical curve shows our approximation 
to the CMS 95% exclusion. The dotted cyan inner (outer) contour shows the 68% (95 %) Bayesian 
credibility region. 



the number of times it appears. We write the point's parameters as 

Pi = {mQ, Aq, mi/2, tan 13, mt,mb{mb),a{Mz),asiMz)}i. 

The Markov chain Monte Carlo statistically gives us proportionality of the weight with the 
posterior probability distribution pciPi) {prior to the inclusion of LHC data) of the point: 
PciPi) = kpG^Pi), where k is some constant which does not depend on pi. 

In order to take into account any new independent data with likelihood Lnew{Pi), 
we have the combined likelihood Lc{pi) = Lnew{Pi)LG{Pi) which combines the global fit 
data with the new data. Thus the combined posterior probability density function is 
Pc{Pi) = cpG{pi)Lnew{Pi), where c is a constant that takes into account the fact that the 
Bayesian evidence changes through the introduction of new data. In terms of the Markov 
chain Monte Carlo point list, we obtain that Pc{Pi) = ckpG{pi)Lnew{Pi)- Ignoring the 
constants, we therefore calculate the new posterior by re-weighting each point, multiplying 
its weight by the likelihood of the CMS ax search in Eq. 2.4. By plotting the posterior 
probability distributions before and after the aT-re-weighting, we can then examine the 
effect of the ax exclusion data on the CMSSM fits. 

We show the fits marginalized over all parameters except for mo and in Fig- 3. In 
order to guide the eye, we have added the 95% exclusion contour (the light, almost vertical 
curve), but we remind the reader that the likelihood was taken into account, not just 
a veto from this curve. Nevertheless, we see in the figure the expected behavior: there is a 
small volume of good fit towards the bottom left of the plot that is "chopped off" by the 
ut exclusion limit. There are also non-trivial effects on the parameter plane: we see from 
Fig. 2a that there is a slight ~ la preference of the data against high mii2, since there was 
an excess of observed events versus expected background. Such effects are present in the 
plot, but are easier to see when the posterior probability density is marginalized down to 
one dimension only (see below). 
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Figure 4. Global CMSSM fits: (a) excluding the CMS ar search and (b) including the CMS 
ut search likelihood. The posterior probability of each bin is shown as the background colour, 
normalized to the maximum bin probability. The dotted cyan inner (outer) contour shows the 68% 
(95 %) Bayesian credibility region. 



We show the marginalization to the ?7ii/2 ~ tan/3 plane in Fig. 4. Here, we start to 
see the effect of non-trivial correlations of the CMSSM parameters from the indirect fits. 
The long, thin vertical region on the left-hand side of the figure corresponds to the ^— pole 
region, where m-^o ~ Tnh/2 and the lightest neutralinos x? efficiently annihilate in the 
early universe through the s— channel higgs {h) pole [42]. This region is at large values of 
mo > 1400 GeV and low mii2 [43]. Examining Fig. 2a, we see that this region is given a 
likelihood enhancement by the ax search, and we see evidence for this in Fig. 4b: it gains 
some 68% level contours as its probability mass increases. Also, the region at high tan /3 
within the 68% level contour has some of its probability mass removed, since it is at large 
values of mi/2> where there is a slight penalty compared to smaller values of ?7ii/2 < 400 
GeV, due to the ut Ic excess of events. 

We display the probability distributions of the individual CMSSM parameters marginal- 
ized over all other parameters in Fig. 5. We see from Fig. 5a that, because of correlations 
between the parameters in the global fit, the ay-search slightly prefers values of mo ~ 100 
GeV. This is because the small mi/2 region is boosted by the la excess in the number of 
observed events, which also has a small mo in order to obtain sufficient neutralino annihi- 
lation in the early universe (in the stau Co-annihilation region). Fig. 5b shows that small 
values of mi/2 are relatively disfavored by the ut search, except for the spike at mi/2 ~ 120 
GeV which corresponds to the /i— pole region and is enhanced by including the ay search. 
Intermediate values of mi/2 of 300 — 500 are also relatively favored compared to larger 
values. The ^o distribution is relatively untouched by the search, as shown in Fig. 5c, 
whereas intermediate values of tan /3 = 8 — 20 are slightly preferred by it: see Fig. 5d. 

Next, we turn to the effect of the ay search on sparticle masses. We plot the posterior 
probability distributions of a representative sample of four different masses in Fig. 6. The 
right handed squark mass m^^ is squeezed somewhat toward intermediate values of 600- 
1200 GeV by the ax search, as shown in Fig. 6a. A similar effect is found in the right-handed 
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Figure 5. Effect of the ar-search on one dimensional probability distributions of CMSSM param- 
eters. The area of each histogram has been normalized to 1 and labeled Tncl. ay' ('Excl. ar') if it 
includes (excludes) the CMS ar search. 



selectron mass rrig^, which is squeezed somewhat toward intermediate values of 160-200 
GeV, as shown in Fig. 6d. m^o and the gluino mass nig are correlated strongly with mi/2 
and Figs. 6b, c show that the overall features follow those of the mi/2 distribution: a spike 
corresponding to the /i— pole at light masses and an enhancement of intermediate values 
of the masses to the detriment of more extreme values. 

The probability distributions of masses found predict a probability distribution for 
the cross-section of sparticle production at the LHC. Approximately 1 fb~^ of integrated 
luminosity is expected to be collected over the next year at the LHC. Thus, a "weather 
forecast" of whether it will be raining sparticles (or whether we shall find a SUSY desert), 
is of special interest. We calculate the total sparticle production cross-section cr suSY 
predicted by our fits at the LHC, with no cuts or modelling of detector effects. This quantity 
then provides a crude upper bound on the actual SUSY signal cross-section that may be 
observed. However, we do not attempt here to quantify the probability of SUSY discovery 
in the next year, since such an inference is not yet robust due to the insufficient constraining 
power of the data. It appears from Fig. 7 that the prospects for SUSY discovery next year 
marginally improve after the inclusion of the ax search results. 



-10- 




Figure 6. Effect of tlie aT-searcli on tlie probability distributions of sparticle masses in the 
CMSSM. Tlie area of each histogram has been normalized to 1 and labeled 'Inch q;t' ('Excl. ar') 
if it includes (excludes) the CMS ar search. 
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Figure 7. Effect of the ar-search on the total SUSY cross-section (tsusy in the CMSSM in pp 
collisions at y/s ~ 7 TeV. The area of each histogram has been normalized to 1 and labeled 'Inch 
q;t' ('Excl. ar') if it includes (excludes) the CMS ar search. 
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4 Summary and Conclusions 



The CMS ot search has significantly extended the previous exclusion limits in the CMSSM. 
We have examined its effect on global fits of the CMSSM to the anomalous magnetic 
moment of the muon, the dark matter relic density and electroweak observables while taking 
into account previous direct searches for supersymmetric particles and Higgs bosons. The 
search nibbles away at the part of the parameter space where both squarks and gluinos are 
light, but also has some other non-trivial effects. Because the search saw a slight excess, 
intermediate masses acquire a small relative preference, with associated effects on tan /3 
because of parameter correlations in the global fit. 

It is remarkable that even with the crudest of detector simulations (consisting of pT and 
r] cuts), and leading-order event generation, our 95% exclusion contour in a certain CMSSM 
parameter plane is close to CMS's more sophisticated treatment. However, in order to be 
sure that we approximately reproduce the CMS likelihood elsewhere in parameter space, 
we require more details from the experimental publication. We strongly advocate that the 
experiments publish more detail in order to allow more detailed checks: at least additional 
confidence level contours or, even better, the likelihood across a parameter plane, or ideally 
a RooStats workspace. 

As this work was being completed, ATLAS produced a search for jets, one lepton and 
missing transverse momentum [44] which extends the exclusion in the CMSSM plane as 
compared to the CMS search that we consider. The ATLAS search has an observed number 
of events lower than the expectation due to SM background (although the measurement 
is statistically compatible with the background). We leave the inclusion of this data to 
a future publication: it is likely to be computationally more difficult to include in our 
fits than the ones contained in the present paper because the signal rate may well have 
a dependence on all of the CMSSM parameters. Event simulation is a CPU-time bottle 
neck, and if it depends on more than just two or three parameters, is likely to need to 
be repeated for different sampled fit points since it couldn't be easily interpolated in a 
reasonable amount of CPU time anyway. The ATLAS data will thus require other more 
computationally intensive techniques to analyze. Some of the effects that we have observed 
from the aT CMS search, such as an enhancement of intermediate values of sparticle masses 
due to the slight excess of events, will presumably be reversed when the recent ATLAS 
data is taken into account. 

A Validation of tan 13 — Aq Independence of the Likelihood 

We now wish to check our approximation that the likelihood C is independent of the 
other CMSSM parameters in the fit aside from rriQ and Mi/2- In order to do this, we 
first take one hundred different sampled points at random from the 2.7 million different 
points in the fits, before the likelihood has been re-weighted by the ay— search. Each point 
has different Standard Model and CMSSM parameters, and the probability of selecting 
it is proportional to its posterior probability density of the fit to indirect data. These 
points therefore represent a faithful sampling from the fit, and we may estimate different 
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Figure 8. Check of the mdependence of the hkehhood with respect to tan/3 and Aq in the CMSSM 
fits, (a) Aq and tan/3 values of the sampled points, (b) binned distribution of Dx^, the difference 
between the Ax^ calculated by our Aq — tan {3 independent approximation and the Ax^ calculated 
by event generation. 



quantities by calculating their probability distributions from this faithful sampling. We 
show, for reference, the values of Aq and tan f3 of the sampled points in Fig. 8a. For each 
of these points, we generate 10 000 SUSY signal events and estimate the numbers of events 
passing cuts in the ax search, as in section 2, turning the number of events into a A^^ for 
the search. We then compare this Ax^ obtained by simulating events with the one from 
the linear interpolation of Aq = 0, tan /5 = 3 in Fig. 2a. The difference between the two, 
Dx^, measures how much our approximation is violated. 

In Fig. 8b, we bin in Dx^, showing how many of the 100 points lie in each bin. Up to 
an overall normalization of 1/100 then, this is an estimate of the probability distribution 
of Dx^ in our fit. All of the 100 points sampled have Dx^ < 0.2 except for six, which are 
shown by asterisks in Fig. 8a. The most egregious has Dx^ = 3.7, but is heavily ruled out 
by the ar-search, with a Ax^ of 15. The others all have niQ < 100 GeV, mi/2 ~ 300 GeV, 
tan f3 £ (5, 10) and Aq £ (—0.5, —1.0) TeV. The fact that the egregious points are clustered 
in parameter space can be understood as follows. The large values of —Aq at small ttiq 
and significantly change the up and down squark masses at these points compared to 
Aq = via the MSSM renormalisation group equations (RGEs)^, having a resultant effect 
on the overall cross-section and hence on the number of SUSY signal events passing cuts. 
Negative Aq is preferred by the global fit as compared to positive Aq because it tends to 
give larger lightest CP-even Higgs masses by enhancing the stop mixing and enhancing the 
(dominant) stop loop correction to the Higgs mass. Once small mi/2 has been selected, the 
global fits preferentially select tan/3 G (5,10), as Fig. 4a shows. The anomalous magnetic 
moment of the muon prefers these values of tan (3 for such light sparticles. For small niQ 
and mi/2, Aq < —1 TeV results in negative stau mass squared parameters, resulting in a 
charge breaking minimum of the scalar potential. Such values are therefore disallowed in 

^The trilinear couplings enter the RGEs for the SUSY breaking masses of the Higgs bosons, which then 
enter the RGEs of the first two families of squarks. 
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the global fit. 

We conclude that the assumed independence to parameters other than mo and 
is quite a good approximation, since, excluding the most egregious point (which is ruled 
out anyway), we obtain Dx^ = 0.04 it 0.22. 
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